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Abstract 

The growth of two-dimensional lattice bond percolation clusters through a cooperative 
Achlioptas-type of process, where the choice of which bond to occupy next depends upon the 
masses of the clusters it connects, is shown to go through an explosive, first-order kinetic phase 
transition with a sharp jump in the mass of the largest cluster as the number of bonds is increased. 
The critical behavior of this growth model is shown to be of a different universality class than 
standard percolation. 
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Percolation concerns the formation of long range connectivity in systems [l], and has 
applications to numerous practical problems including conductivity in composite materials, 
flow through porous media, and polymerization [2]. In the usual random (Bernouilli) per- 
colation, bonds are formed randomly and independently throughout the system, and at a 
critical concentration of occupied bonds, a finite fraction of the sites (vertices) are con- 
nected together and percolation takes place through a continuous or second-order transition. 
By the universality hypothesis, all large-scale properties near the transition point (i.e., the 
fractal dimension, critical exponents) are the same in a given dimensionality, irrespective of 
the microscopic details of the rnodel 

Recently, Achlioptas et al. [3] modified the growth of percolation clusters to produce 
a first-order kinetic transition on the mean-field-like random graph, through a procedure 
that is known as an Achlioptas process. In this optimization process, introduced to study 
problems in graph theory, two alternate choices of adding a bond are considered, and a 
specific strategy to favor a desired result (such as delaying or accelerating the appearance of 
a giant component, or the formation of a Hamiltonian cycle) is followed QQ. InRef. 0, the 
authors considered a process of growth of clusters in the Erdos-Renyi random graph network 
model in which two unoccupied edges between clusters are chosen at random, and the one 
that minimizes the product of the two connecting cluster masses is preferentially chosen as 
the next occupied bond. This process led to a sharp jump in the growth dynamics, which 
the authors labeled as explosive percolative growth, in contrast to the normal percolation 
phase transition, on both regular and random lattices, where the transition is second-order 
and continuous. 

Erdos-Renyi random networks are formed by randomly linking pairs of points, irrespective 
of their distance apart. These networks are essentially mean-field and infinite-dimensional, 
and for this problem, only the mass of the clusters needs to be kept track of. There is a 
close connection between the percolation on Erdos-Renyi networks and percolation on the 
infinite-dimensional Bethe lattice and polymerization ^f non-looping branched polymers, 
whose percolation or gelation theory goes back to Flory [6] . When bonds are added randomly 
between sites in the random graph, the net effective probability that a cluster of mass Si 
and a cluster of mass S2 are joined is proportional to the product of their masses S1S2. This 
leads to the product kernel of the polymer growth process as described by the Smoluchowski 
equation, and this growth process is equivalent to percolation on the Bethe lattice 
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Choosing the bond with the minimum product leads to the minority product rule (PR) of 
Ref. y. 

While random and scale-free networks have recently received a great deal of attention 
sl, many actual networks of interaction, whether physical or social, are restricted spatially 
and bonds can only form between close neighbors. The two-dimensional percolation model 
is the strongest form of such a restriction and has been the subject in intense study for over 
fifty years p|. Random percolation can be viewed as the result of a growth process where 
bonds are added to a single cluster (as in the Leath algorithm) or between different clusters, 
as considered in [9|]. Both of these processes lead to second-order transitions, as do their 
more dynamic cousins, the directed-percolation and contact processes. Here we consider the 
question of whether the PR rule produces a discontinuous transition in the two-dimensional 
percolation model also. Indeed, we find that the transition is very sharp and explosive, and 
find strong indication that the transition is indeed first-order. 

We consider bond percolation on Lx L square lattices with periodic boundary conditions 
in both directions, with n = L?' sites. Bonds are added randomly and one at a time, keeping 
track of the current cluster structure of the system (the so-called Newman-Ziff process) jol, 
modified by the consideration of two unoccupied bonds that could connect distinct clusters. 
The bond that minimizes the product of the masses of the two clusters it joins is preferentially 
chosen to become occupied, and the clusters are merged into one. Time t represents the 
number of successful bonds added, and each new bond, which always connects different 
clusters, reduces the number of clusters by one. Sites with no bonds are considered to be 
clusters of one site, so initially we have n clusters. After t bonds are added, the number of 
clusters N in the system equals n — t. We characterize the mass of a cluster by the number of 
sites, s. Because we only add bonds between sites of different clusters, the bonds themselves 
form "minimum spanning trees" over each cluster, yet for the conventional growth process, 
where bonds are added one at a time, standard percolation clusters (characterized by the 
sites that are connected ) are created. 

We have carried out simulations of the cluster growth using both the regular and PR 
processes, like in Ref. ^| measuring the maximum cluster size C as a function of t. The main 
plot of Fig. [1] shows the results of these simulations for single realizations on a 1024 x 1024 
lattice. In the PR model, the transition is quite sharp. An expansion of this curve shows that 
the jump substantially occurs over a small number of time steps, as shown in the inset for 
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FIG. 1: (color online) Plot of regular percolation (blue) and the PR (red) for bond percolation on 
a lattice of size 1024 x 1024, showing the delayed and explosive growth in the PR model. Points are 
every 1024 time steps. Inset shows PR model on an expanded scale, on systems of sizes L = 256, 
1024, 2048 and 8192, with the scaled axes A/n vs. t/n with n = L"^. 

different system sizes (here with the axes scaled by n). In contrast, for regular percolation, 
the transition is much smoother. The simulations on different size lattices show that the 
first-order transition in the PR model is robust. 

For regular percolation, we can locate the transition point exactly. While we don't know 
the effective overall bond occupation fraction p, which includes bonds between sites in a 
cluster beyond the spanning tree, and thus don't know when the square-lattice threshold 
Pc = 1/2 is reached, we can identify the transition by finding the point where the number 



of clusters per site reaches its critical value lOl. Ill] 



Nc 3^3 - 5 



n 



0.098076 , 



(1) 



which is known to have small finite-size corrections for a periodic system llj. This density 
will be reached when t/n = l- Njn = (7 - 373)72 ^ 0.901924. Of course, for a normal 
system at pc, there will normally be some fluctuation in Nc-, but for large systems these 
fluctuations will be small, so finding where the density exactly equals 1 — N^/n gives an 
excellent estimate of the critical point. We have verified that, at this point, the average C 
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FIG. 2: (color online) Scaling of A/n = (to — ti)/n for regular percolation (upper curve) and the 
PR model (lower curve) with system size L. 

for different size systems scales as with D ^ 1.8953, consistent with the theoretical value 
91/48 ~ 1.8958. On the other hand, for the PR rule, we have no a priori knowledge on 
where the transition point should be. 

Achlioptas et al. characterized the transitions by two times: the time Iq where C equals 
a/ti, and the time ti where C equals n/2, and then considered A = ti — t^. For unbiased 
growth, they argued that A should be proportional to n, while for the PR growth they found 
A ^ n^f^ . We have measured the same quantities for regular the PR growth on square 
lattices with L = 32, 64, . . . 8192, with the number of realizations ranging from 100000 for 
the smaller sizes to 100 for the largest. The results for A are plotted in Fig. [2l 

This figure provides clear evidence that the two transitions are of a quite different nature. 
The regular percolation model does not correspond to A/n going to a constant as in the 
Erdos-Renyi case, but instead decreases as A/n ~ /^-o.sss _ ^^-0.192^ fact, ti converges 
rapidly to tc (because, for finite lattices of this size, at the critical point C/n is close to 
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1/2, which is the condition used to calculate ti). The main variation in A is due to the 
variation in to- By normal scaling arguments, we expect the typical cluster mass s* to scale 
as s* ~ Ip — Pcl^^^'^ with a = 36/91 ~ 0.3956 in 2d l|], so setting s* = ^/n = L to correspond 
to the point to, and using the fact that po — pc is proportional to (to — tc)/n ^ —A/n for 
small Po — Pc, we deduce that 

A/n~n-'^/2 = L'=^^/9i , (2) 

which agrees fairly well with our numerical observations. Actually, the maximum cluster 
size C differs from the typical cluster size s* by terms logarithmic in L, which would add 
logarithmic corrections to the above prediction and may account for the small deviations 
seen in the behavior from ([2]). 

We emphasize that po — pc is proportional to, but not equal to, (to — tc)/n, because in 
our growth process, t = n ~ N reflects the number of bonds linking different clusters added 
to the system, and does not include bonds within existing clusters. That is, p is not equal 
to t/n. The number of clusters N{p) = n — t has a high-order singularity \p — pd"^^"" with 
a = —2/3 at Pc, but this is masked by the lower-order continuous terms. Thus, it is valid 
to assume (to — tc)/n oc po — Pc to first order. 

For the PR model, we find the much different behavior: A ~ /^-o.esa Closer examination 
shows that both to — tc and ti — tc scale with the same exponent of about —0.70, with 
coefficients —1.429 and 0.040, respectively. Thus ti/n changes very little with L, and has 
a well-defined limiting value ~ 0.9925 as L — oo. This limiting value implies a critical 
density of clusters Nc/n ^ 0.0075 at the transition point, much lower than the value Nc/n 
for regular percolation given in Eq. ([T]). 

We illustrate the critical behavior of the two models in Fig. [3l where we show snapshots 
of the surface for each model for C ^ 72/2 on a 128 x 128 lattice. Here we show only the 
sites of the lattices, colored differently according to the clusters they belong to. The lower 
number of clusters in the PR model is clearly evident. We also show in Fig. |l]a closeup of 
the PR case, showing the connecting spanning-tree bonds. 

This marked difference in the critical behavior between the PR model and regular perco- 
lation implies an underlying difference in the nature of the transition. The behavior in Fig. [1] 
points to that transition being of first order in the sense that once the large clusters form, the 
process rapidly goes to near-completion and the maximum cluster size has a discontinuous 
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FIG. 3: (color online) Regular (left) and PR (right) percolation models at C ~ n/2. 




FIG. 4: (color online) PR model close-up, showing the spanning-tree bond structure. 



nature. 

Having a simple lattice growth model showing strong first-order behavior is useful in 
modeling variety of discontinuous and explosive growth processes. Kinetic growth processes 
have been the object of extensive study over the last two decades yet this type of growth 
behavior is unusual in simple models. An example of a growth process that leads to behavior 
almost identical to the PR curves shown in Fig. [1] is in the dynamics of a catalyst system 
(oxidation of carbon dioxide) with defects randomly scattered on the catalyst surface 13|. 
The defects have the effect of creating regions of different sizes, and as the CO partial pressure 
increases, regions "poisoned" with CO coalesce along paths of least resistance and create a 
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very similar explosive growth of surface coverage. Many questions for the present model, 
such as fractal properties of the resulting clusters, other critical exponents, dependence upon 
growth rules, behavior in higher dimensions, etc., remain. 

The author thanks Dietrich Stauffer for encouragement and many helpful comments, 
and acknowledges support from the U. S. National Science Foundation grant number DMS- 
0553487. 
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